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With an eye towards applications in chemistry, we characterize Wigner quasiprob- 
abihties for boson-type and fermion-type wave functions or density matrices in turn. 
Our characterization appears to be more convenient that the one given by O'Connell 
and Wigner [1] some time ago. 

1 Introduction 

For eighty years, Wigner quasiprobability distributions (or Wigner functions for short) have 
been an important tool of quantum physics, successfully adopted in statistical mechanics [2], 
optics [3] and now quantum chemistry [4]. 

By definition, a pure n-body Wigner quasiprobability distribution in terms of the wave 
function \E' is given by 



with the notation ) for n bodies, where Xi G M.'^ (say, d = 3 for ordinary 

space), and similarly for p and z. Here \l/ denotes the wave function on momentum space 
corresponding to We have taken units so that h = 1. 
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(1) 



In this context, quasiprobability means that averages of this expression with classical 

phase-space observables reproduce the expected values predicted by standard quantum me- 
chanics. Wigner functions can describe statistical mixtures as well, by the usual convex 
linear combinations. 

It is not easy to characterize Wigner functions, although necessary and sufficient condi- 
tions (not involving wave functions explicitly) for a phase-space function to be an admissible 
Wigner quasiprobability distribution are known [4-6] . 

Now, a natural question is to determine when a Wigner quasiprobability distribution 
does correspond to a wave function symmetric or antisymmetric under permutations of its 
variables. This was posed since the early days, mainly in view of applications to kinetic 
theory and statistical physics: see the references in [1]. But only the latter article purported 
to offer a general answer.^ 

On their characterization, O'Connell and Wigner do write: "It must be admitted that 
this equation, for the distribution function, postulating the Bose statistics for a system of 
spin particles, is much more complicated that the corresponding equation for the density 
matrix". Et pour cause: formulas (14a) or (14b) in [1], together with formula (11) look 
unwieldy, so much as perhaps to have discouraged borrowing of the more interesting part 
of their paper [1, Sect. IV], on systems of spin 1 particles, which unfortunately we have not 
seen exploited elsewhere. No examples are given. 

The development of density functional theory based on Wigner quasiprobabilities [4] 
obliges us to revisit the matter here, from the standpoint of quantum chemistry. This 
means that our particles are electrons, that is, fermions; however, the orbital parts of the 
corresponding wave functions may be symmetric or antisymmetric. For closed-shell config- 
urations, the spin-free part of the ground state is symmetric, while that of the first excited 
state is typically antisymmetric. 

It should be clear that, whether ^ is symmetric or antisymmetric, the function (1) is 
symmetric under particle interchange, 

W^{Xi, . . . ,Xi, . . . ,Xj,. . . ,Xn;Pl. ■ ■ ■ ,P,. ■ ■ ■ -Pj,. ■ ■ ,Pn) 

— W\^[Xi, . . . , Xj, . . . , Xi, . . . , Xn] Pi, . . . , Pj, . . . , Pi, . . . , Pn)- 

We readily give a long-due simple answer to the physical indistinguishability question for Wi;,. 
Our characterization takes the form of a mere preservation/change of sign under (the same) 
permutation of two variables, respectively for the symmetric/antisymmetric case — ^just like 
in the ordinary formalism of quantum mechanics. This makes it trivial that the square of 
a permutation induces the identity, which is not so obvious in [1]. A last bonus is that our 
conditions are akin to those of the general characterization problem for Wigner functions, 
as expounded in [4, App. B]. We illustrate them with a few example classes. 

Consider first 2-body functions. Bringing in mean and difference coordinates, or, in 

^We put aside the second- quantized approach to Wigner quasiprobability, which has drawn scant atten- 
tion. 
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chemists' jargon, extracule and intracule coordinates, respectively given by 

^{xi + X2), r = ^{xi-X2), (2) 

with the customary abuse of notation, the symmetry/ antisymmetry conditions (say, on con- 
figuration space) ^{xi, X2) — ±^{x2, Xi) read: 

It is not hard to see that with 

P = ^iPl+P2), P=-^(P1-P2), (3) 

the meaning of W{R, r; P, p) is unambiguous. This is due to the hnear symplectic invariance 
of the Wigner function formahsm. 

If instead of the wave function we employ the density matrix p{xi ), the same 

transformations give rise to 

p{R, r, R, r ) = ±p{R, ±r, R, r') = ±p(i^, r, R', ±r'). 

In this case, Wp{x; p) = tt^'^" J p(x — z; x + z) e^**''^ dz. 

Since the discussion will turn around the intracule variables, it is worth regarding R, P 
as parameters and introducing the following notation: 

c^H,p(r";p):=Pr(i2,r;P,p). 

Let us denote the following partial Fourier transform on the intracule set of variables: 

ujR^p{v;p) := j ujR^p{r;p) e^'""^ dr = ujr,p{-v; -p). 

The last equality holds if ^ or p possesses the symmetry or the antisymmetry property. 

2 When does W come from a symmetric wave function? 

Proposition 1. A Wigner two -body function comes from a boson-type density matrix if and 
only if, for all v and p: 

^^rA'"'^p) ^^rAp'^'")- (4) 
Proof. Consider first the following integral: 

CjRp(v;p)^Ar; I ^{R- Z,r - z)^*iR+ Z,r + z)e^'^-^+'^'P-^ e^'"'-'' drdZ dz 

^2d J 

= ^ J ^{R- Z,z-r)^*{R+ Z,z + r) e^^^-^+^i'"^ e^'P'^ dr dZ dz 



I- 
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Conversely, retracing our steps till the second equality, one verifies that the identity (4) holds 
only if ® ^* defined on x M.'^ is symmetric. 

For the general case of statistical mixtures, exchanging the second (say) set of variables 
we find: 

ojRAf^-, P) = ^ J p{R- Z,r-z;R+Z,r + z) e^^^'^+^^P'^ e'^^"" dr dZ dz 
= ^ J p{R- Z,z-r;R+ Z,z + r) e^'^-^+^'''-'- e^'^'^ dr dZ dz 

As before, this will hold only if p is symmetric. We could have argued as well from convex 
linear combinations of pure Wigner functions. □ 

The reader may check, after a few steps, that (4) is mathematically equivalent to (14b) 
together with (11) in [1], when d — 3. 

In quantum chemistry based on the Wigner function formalism [4] , use of Gaussian basis 
sets is if anything more natural than in standard quantum chemistry [7] . This motivates the 
example that follows. 

Example 1. Take as a boson- type wave function the symmetric product of two general 
Gaussians centered at the origin: 

where, for j = 1,2: 

V'j(^) = ^ e"^'^^-^'-^^^'^^-^" with dj > 0, bj e R. 
(The normalization factor C fulfils 



Cr^ - 2 + 2. 



1^/1^/2 



(dl + d2)^ + {hdi - h2d2Y ' 

but this is unimportant here.) The corresponding two-body quasidensity is 

iy(Xi,X2;Pl,P2) OC Wii{Xi;Pi)W22{x2\P2) + W22{Xl\ Pl)Wii{x2\ P2) 

+ Wi2{Xi]Pi)W2l{x2]P2) + W2l{Xi]Pi)Wi2{x2]P2)- (6) 

Here Wj^ represents an interference, namely, 

^1/4 >l/4 

djk ■■= ^{dj + dk) + ^{bjdj - bkdk), bjk ■= \{bjdj + 6fc4) - '^{dj - dk), 
Ajk := djk + bjf./ djk, Bjk := bjk/ djk- 



4 



The quadratic form in the exponent of the Gaussian is given by a symmetric, symplectic 
matrix with positive definite real part. When k — j, we have a Gaussian pure state, 



whose coefficient matrix is real, symplectic and positive definite [8]. 

To see that the quasidensity (6) fulfils (4), change variables according to (2) and (3), and 

let 



exp{-Ajkxl - 2BjkXiPi - dj^pf) exp{-AkjX2 - 2BkjX2P2 - d^jPi 



kdkj 

' d 

^ exp [- Re Ajk {R^ + r^) - 2i Im AjkRr - Re dj^ (P^ + p^) - 2i Im dJ^Pp] 



TT^ y djkdkj 
1 / didk 



T^^Md^kdkj ^ ; ^.x.^.-,.^.- -^^j ^.x.^.-,-, 

X exp [-2 Re Bjk{RP + rp) - 2i Im Bj^irP + Rp)] . 
Now, multiplying by e^*^*" and integrating with respect to r, we obtain 
V27r-^/^ 



exp [- Re AjkR^ - Re dj^ (P^ + p^) - 2i Im dJ^Pp] 



^ dj + b'jdj + dk + bldk 

i Im AjkR + Re BjkP + i Im Bj^P — iv) 



X exp [—2 Re Bj^RP — 2i Im Bj^Rp] exp 
We have used here that 



Re^ 



jk 



djkdkj{Ajk + Akj) = djdk{dj + b'^jdj + 4 + bldk). (7) 

Now let Wjj{xi; pi)Wkk{x2', P2) —'■ ^jjkk{R, t^', P,p)- With a little more work, using (7) again, 
one checks that 

J \jMR^r;P,p)e''''''dr^ J Xjkkj{R:r;P,v)e'''P'-dr with A; ^ j, 
thereby verifying condition (4) for this example. 

3 When does W come from an antisymmetric wave function? 

Proposition 2. A Wigner two-body function comes from a fermion-type density matrix if 
and only if for all v and p: 

oj R,p{v;p) = -ujR,p{p;v). (8) 
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Proof. We argue from wave functions. As in the boson-type case, we consider the following 
integral: 

URp(v;p) = ^ [ ^{R- Z,r- z)^>\R+ Z,r + z)e^'^-^+^'P-^e^''' '' drdZdz 

J 

= J <i/{R- Z,z-r)^*{R+Z,z + r) e2^^-^+2^^-^ e'^^"^ dr dZ dz 
= - y" ^R,p{z] v) e^'P'^ dz = -Cor^p{p] v). 

Therefore (8) holds if and only if \E' defined on x is antisymmetric. □ 

As a consequence of formula (8), a necessary condition for a Wigner function to come 
from an antisymmetric wave function is u}iip{p,p) = 0. This is not sufficient, since any 
one-dimensional wave function (say) of the form ^'(i?, r) = f{R)g{r) with J g{r) dr = will 
give rise to a Wigner function having that property. 

As for examples, it is clear that Gaussian sets like the ones in (5) with a minus instead of 
a plus sign will do. However, it is more instructive now to take a different tack. In the early 
years of Quantum Mechanics, as a prolegomenon to calculating the energy levels for helium, 
Hciscnberg [9] studied the harmonium, an exactly integrable analogue of a two-electron atom. 
It exhibits two fermions interacting with an external harmonic potential and repelling each 
other by a Hooke-type force. Being simple, but not trivial, this system has been borrowed in 
many contexts. It is sometimes called the "Moshinsky atom" , since Moshinsky reintroduced 
it with the purpose of studying correlation energy [10,11]. Also, it has been employed in 
approaches to Bose-Einstcin condensation [12] and black-hole entropy [13]. 

The harmonium Hamiltonian in Hartrcc-likc units is given by: 

H{Xi, X2]Pl,P2) = ^ + ^ + -{\Xi\' + \x2r) X2\'. 

Introducing extracule and intracule coordinates and the frequencies u := y/k and /i :— 
\Jk — S, the Hamiltonian is rewritten as that of two independent oscillators: 

^ ^ ^ |PP U^\R\^ |p|2 

H ^ Hr + Hr -.^ - — - H — - + — + . 

^ 2 2 2 2 

Example 2. Since the harmonium problem factorizes completely, we may work in dimen- 
sion one. The orbital part of such an eigenfunction is written 4>n{R)ipm{'>^) , where the par- 
ity of ijjmir) must be even for spin singlet states and odd for triplet states. The Wigner 
quasiprobabilities associated to those eigenvectors have the general form: Wn{R, P)Wm{r,p), 
where, with L„ denoting the n-th Laguerre polynomial [14, 15]: 

TT 



TT 
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The generating function of the Laguerre polynomials is: 

oo 
m=0 

Defining 

we obtain 

oo „ 

^-^ 7r(l — x) J 



7r(l — x) 



Thus, r^(t',p) = —Tm{p,v) if m is odd and r^(u,p) = rm(p, i') if m is even. In particular, 
whenever the function = (pnipm is antisymmetric, the corresponding Wigner function 
Wnm{R: P',f:P) — Wn{R, P)Wjn{r, p) in agreement with (8) must satisfy: 

J Wnm{R,P;r,p)e^''''dr = - J Wnm{R, P;r,v) e^'P'dr. 

Needless to say, we could have formulated our conclusions in a completely analogous 
manner in terms of ujR,p{r, s) := J UR^p{r;p) e~^'®'^dp instead of ur^p{v]p). 



4 The general n-body case 

Our results here easily extend to ra-body systems just by considering intracule and extracule 
coordinates for each pair of particles. This gives rise to n{n — l)/2 conditions like the ones 
above. The details are obvious, and we omit them. 



Acknowledgments 

We are most grateful to Joseph C. Varilly for a careful reading of the manuscript. CLBR has 
been supported by a Banco Santander scholarship. JMGB has been supported by a grant 
from the regional government of Aragon. He owes as well to the Zentrum fiir interdisziplinare 
Forschung, for support and warm hospitality. 



7 



References 

[1] R. F. O'Connell and E. P. Wigner, Phys. Rev. A 30 2613 (1984). 

[2] E. P. Wigner, Phys. Rev. 40 749 (1932). 

[3] U. Leonhardt, Essential Quantum Optics, Cambridge University Press, Cambridge, 2010. 

[4] Ph. Blanchard, J. M. Gracia-Bondi'a and J. C. Varilly, Int. J. Quant. Chem. 112 1134 (2012). 

[5] F. J. Narcowich and R. F. O'Connell, Phys. Rev. A 34 1 (1986). 

[6] J. M. Gracia-Bondia and J. C. Varilly, Phys. Lett. A 128 20 (1988). 

[7] F. Jensen, Introduction to Computational Chemistry, Wiley, Chichester, 2007. 

[8] R. G. Littlejohn, Phys. Rep. 138 193 (1986). 

[9] W. Heisenberg, Z. Physik 38 411 (1926). 

[10] M. Moshinsky, Am. J. Phys. 36 52 (1968). 

[11] K. Ebrahimi-Fard and J. M. Gracia-Bondi'a, J. Math. Chem. 50 440 (2012). 

[12] R. Seiringer, Ph. D. thesis, Vienna, 2000. 

[13] M. Srednicki, Phys. Rev. Lett. 71 666 (1993). 

[14] C. L. Benavides-Riveros, J. M. Gracia-Bondia and J. C. Varilly, "The lowest excited configuration of 
harmonium" , physics. chem-ph/1205. 2038. 

[15] C. L. Benavides-Riveros and J. C. Varilly, "Testing one-body density functionals on a solvable model", 
physics.chem-ph/1207.3897. 



8 



